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In this work, the transition between diffusion- limited and ballistic aggregation models was revisited 
using a model in which biased random walks simulate the particle trajectories. The bias is controlled 
by a parameter A, which assumes the value A = (1) for ballistic (diffusion-limited) aggregation 
model. Patterns growing from a single seed were considered. In order to simulate large clusters, 
a new efficient algorithm was developed. For A / 0, the patterns are fractal on the small length 
scales, but homogeneous on the large ones. We evaluated the mean density of particles p in the 
region defined by a circle of radius r centered at the initial seed. As a function of r, p reaches the 
asymptotic value po(A) following a power law p = po + Ar~"< with a universal exponent 7 = 0.46(2), 
independent of A. The asymptotic value has the behavior po ~ |1 — Al' 3 , where (5 = 0.26(1). The 
characteristic crossover length that determines the transition from DLA- to BA-like scaling regimes 
is given by £ ~ |1 — \\~" , where v = 0.61(1), while the cluster mass at the crossover follows a power 
law ~ |1 — A|~ a , where a = 0.97(2). We deduce the scaling relations j3 = 1/7 and /3 = 2v — a 
between these exponents. 
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I. INTRODUCTION 

The pattern formation in noncquilibrium processes is 
a longstanding problem in Statistical Physics 0,0,13 In 
special, the diffusion-limited aggregation model (DLA) 
[4j is a noteworthy example in which a very simple algo- 
rithm generates disorderly fractal clusters. This model 
was related to several physical and biological applica- 
tions, such as electrodeposition viscous fingering 
bacterial colonies 0, and neurite formation U I n DLA 
model, particles released at a point distant from the clus- 
ter execute random walks until they find a neighbor site 
of the cluster and irreversible stick in this site. If the ran- 
dom walks are replaced by ballistic trajectories at ran- 
dom directions, we have the ballistic aggregation (BA) 
model In contrast to DLA, the BA model gener- 
ates disordered nonfractal clusters with nontrivial scaling 
properties @,E3- 

Due to its importance as a fundamental model.sev- 
eral generalizations of DLA model were proposed P, • 
In particular, those models in which the particle tra- 
jectories are biased random walks were investigated 
[H 13 13 III 13 HD- In these models, on the short 
length scales the particle trajectories are common ran- 
dom walks with fractal dimension 2.0, whereas on the 
longer length scales the bias becomes dominant and the 
dimensionality of the walk is 1.0. Clusters grown us- 
ing this type of walk must behave like DLA model on 
the short length scales, while nonfractal patterns are ob- 
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served on the longer ones. Consequently, the mass of a 
cluster of size I is given by 

M(l) = l d tf(l/& (1) 

in which 

, . J const, if x <C 1 / x 

■ t[x) ^ \x d - d f ifx^>l. {Z) 

Here, d is the space dimension, df is the DLA fractal 
dimension, and £ is the crossover radius from DLA- to 
BA-like scaling regimes. This idea was firstly consid- 
ered by Meakin pd| . In his model, the simulations start 
with a single seed at the center of a square lattice and 
the drift of all trajectories is in a fixed lattice direction. 
Along the walk, the particle is moved one lattice unity 
in the drift direction with probability P, or moves to one 
of its four next-neighbor sites with probability 1 — P. 
The model generates patterns with a growth tendency 
in the opposite direction of the drift. The author ar- 
gues that the crossover from the DLA-like structure on 
the short length scales to a compact structure on the 
longer ones is characterized by a length £ ~ P^ 1 . How- 
ever, using a renormalization group approach, Nagatani 
found £ - p-i/(d-d f ) Kim et d [i^lil studied 

lattice models with a global drift to the seed, in which 
the particles have a higher probability to move to the 
nearest neighbor representing the shortest distance away 
from the seed. The pattern morphologies are ruled by 
the lattice anisotropy and their fractal dimension is 1. 
Nagatani 0] considered the effects of positive and neg- 
ative radial drifts in DLA model. In the positive case 
the cluster fractal dimension is asymptotically 2, while 
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eccentric patterns with dimension 1 were found for the 
negative case. Other models 0, 0] consist of the de- 
position processes on a d-dimensional substrate in which 
the walk drift is the substrate direction. 

In the present paper, we are interested in the transi- 
tion from DLA to BA models when the random or ballis- 
tic trajectories of DLA and BA models, respectively, are 
replaced by biased random walks with a random drift di- 
rection. The central concern of this work is the fact that 
all real fractals exhibit scaling only on limited ranges and, 
consequently, a quantitative analysis of both experiments 
and simulations demands a deep understanding of these 
crossovers. The outline of the paper is the following. In 
Sec. [H] the model and the respective computational algo- 
rithm are described. In Sec. IIHI the simulational results 
are presented and discussed. Finally, some conclusions 
are drawn in Sec. IIVI 



II. MODELS AND METHODS 

As in the original DLA and BA models, at the be- 
ginning of the simulations a unique seed localized in the 
center of the lattice constitutes the cluster. Then, parti- 
cles are sequentially released at a circle distant from the 
cluster and execute biased random walks. The distance 
between the center of the lattice and the launching circle 
is denoted by i?/. The biased walks are defined by 

x n+ i = x n + cos(<^ + \0 n ) 

Vn+i = y n + sin(v + A0„), (3) 

where x n and y n are the particle coordinates at the nth 
step of the walk, <p is a random angle that defines the 
bias direction, A € [0, 1] is the parameter that controls 
the random component of the trajectories, and 9 n is a 
random direction. The variables <p and 8 n are in the 
range [— 7r,7r]. Notice that <p is defined at the beginning 
of the walks, whereas 9 n assumes random values for each 
walk step. One can see that the particular cases A = 
and A = 1 recover the BA and DLA models, respec- 
tively. If the particle visits a site neighboring the cluster 
it irreversibly joins to this site. However, if the distance 
between the particle and the cluster is too large, i. e., 
larger than a killing radius Rk, the particle is excluded 
and a new one is released at the launching circle. In order 
to determine when a walker is neighboring a cluster site, 
its lattice position was defined as the next integer value 
of its real coordinates defined by Eq. 0. The majority 
of the results presented in this work refer to the square 
lattice version of the model. 

The previously introduced variables (Ri and Rk) must 
be as large as possible. However, computational limita- 
tions restrict the use of Ri and Rk values. We defined 
Ri = Rmax + Rq, where R ma x is the maximum distance 
from the center of the lattice of a particle belonging to 
the cluster. For DLA and BA models R can be of the 
order of some lattice units However, for the model 




FIG. 1: Schematic representation of the model. 



with biased random walks the patterns morphologies are 
strongly dependent on this value. Our tests suggest that 
the patterns become insensitive to Rq variations when 
Ro > 300, in agreement with the values adopted by Kim 
|16j for a model of deposition of biased random walks 
on a substrate. Thus, R = 400 were used in all simu- 
lations. The killing radius Rk must be 10 to 100 times 
Rmax for very large DLA clusters Q whereas a Rk only 
some lattice units larger than R ma x is necessary for the 
BA model. Due to the bias present in the random walks, 
we used the same strategy adopted by Kim 0], i. e., 
Rk = 2R max + Rq. Fig. ^ illustrates two tentatives, one 
successful and the other frustrated, to add a new particle 
to the cluster. 

To analyze the transition between BA and DLA it is 
necessary to simulate large clusters using lattices con- 
taining up to 10 4 x 10 4 sites, specially when A < 1. Con- 
sequently, the computational time becomes prohibitive 
and an efficient algorithm is necessary. A technique com- 
monly used to simulate large DLA and related models is 
to allow the particles execute long steps at random direc- 
tions if they are far from the cluster [l8l Il9l I20L l2*lj . This 
procedure is correct because the probability of a random 
walker crosses the circle centered on its initial position 
in a given angle <j> is uniformly distributed in the interval 
[— 7r,7r]. However, for the biased random walkers this is 
not true. Indeed, the probability density distributions 
are concentrated around the direction (p. For A not very 
close to the unity or for large steps, the probability distri- 
butions are very well fitted by Gaussian curves centered 
at cp, as illustrated in Fig. [5] Using this fact, the fol- 
lowing procedure was adopted. If the distance between 
the random walker and the cluster is larger than a value 
R s + 5, it executes a jump of length R s . Thus, a long 
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FIG. 2: Angle distributions for biased random walkers with 
preferential direction tp — and the correspondent Gaus- 
sian fits. The curves correspond to A = 0.99 and R s = 200 
(squares), A = 0.99 and R s = 1000 (triangles), and A = 0.90 
and R 3 = 200 (circles). Other parameter sets used in the sim- 
ulations provide fits better than those shown in this figure. 
For all couples A and R s used in the simulations a correlation 
coefficient r 2 > 0.999 was required. 



jump can not lead the walker to a distance smaller than 



(5. The jump direction is (j) — <P ~ 



dom number between 
distribution 

P(4>g) = - 



where <b„ is a ran- 



-7r and 7r selected from a Gaussian 



: exp 




(4) 



In order to obtain the Gaussian width a for each cou- 
ple A and i? s , a large number of biased random walks 
(10 6 — 10 7 ) were simulated and a histogram of the prob- 
abilities built (Fig. |2| . Then, the a value can be deter- 
mined using least square Gaussian fits. The quality of 
fits is improved as larger R s values and smaller A values 
are used. Therefore, the R s values should be sufficiently 
large for reproduce good fits, in special for A < 1. Fur- 
thermore, several R s values can be used in a same sim- 
ulation improving the algorithm efficacy. We used two 
values: R s = 200 and R s = 1000. The a values used in 
the simulations are shown in table Q] Also, 5 = 20 were 
used. All tests show that the growth patterns are not sen- 
sible to the 5 value. It is worth to note that the Gaussian 
distribution is not normalized in the interval [—it, tt] and, 
obviously, this is not the actual angle distribution for the 
present problem. However, the very good fits to the angle 
distributions justify the use of the Gaussian functions. 



III. RESULTS AND DISCUSSIONS 

The first stage of the present work was to confirm the 
validity of the previous defined algorithm. We simulate 



0.100 
0.300 
0.500 
0.700 
0.900 
0.950 
0.990 
0.995 



0.02585 
0.07591 
0.12535 
0.18206 
0.31933 
0.45109 
1.03575 
1.53278 



0.01202 
0.03413 
0.05613 
0.08139 
0.14216 
0.20057 
0.45130 
0.63451 



TABLE I: a values determined with the Gaussian least square 
fits. 



relatively small lattices containing 10 3 x 10 3 sites with 
and without the optimization for A = 0.99. In Fig. [3] 
comparisons between clusters generated with (top) and 
without (bottom) the optimization are shown. Compar- 
ing the patterns, one can see that they are statistically in- 
distinguishable. Using the mass-ratio method, the fractal 
dimensions of the patterns generated with and without 
the algorithm were d f = 1.70±0.02 and dj = 1.72±0.02, 
respectively and the exponents of the radius of gyration, 
defined by r g ~ rfi (n is the number of cluster parti- 
cles) 0, were < = 0.560 ± 0.003 and C = 0.561 ± 0.002, 
respectively. These exponents reinforce the algorithm va- 
lidity. Concerning the computational time, a single run 
to generate one of the clusters shown in Fig. [3J without 
the optimization takes about 1 hour in a 3 GHz Pentium 
IV, but the same simulation is done in 10 minutes us- 
ing the optimization. Therefore, even for small lattices 
the simulation performance is greatly improved when our 
optimized algorithm is used. 




FIG. 3: Two clusters generated (a) using or (b) not the opti- 
mization. In these simulations, lattices of size L = 1000 and 
A = 0.99 were used. 
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Figure 21 shows growth patterns for distinct A values. 
These patterns were generated with the optimization 
at lattices containing 10 3 x 10 3 sites. The simulations 
stopped when the aggregate reaches the lattice edge. A 
continuous transition from disordered and dense to ram- 
ified clusters is observed. For small A values the pat- 
terns are essentially BA-like but, the patterns become 
very similar to the DLA clusters as A — ► 1. Indeed, the 
cluster generated with A — 0.995 is characterized by the 
square lattice anisotropy, a signature of the DLA model 
PI Q. However, one expects that all patterns become 
asymptotically homogeneous with a finite characteristic 
size for the empty regions. 

In order to quantify the DLA to BA morphology tran- 
sition, the mean particle density in the inner regions of 
the cluster was evaluated. This mean density p(r) is de- 
fined as the ratio between the number of occupied sites 
and the total number of sites in a region delimited by a 
circle of radius r centered at the initial seed. Since one 
expects asymptotically nonfractal clusters, the density 
must reach a finite value po as r — > oo. Nevertheless, the 
approach to the constant density is very slow and takes 
a scale invariant form 



p(r) = po + Ar 1 . 



(5) 



Here, 7 is a correction to the fractal dimension and A a 
constant. This scaling hypothesis was also used by Liang 
and Kadanoff 10] to study the driven ballistic aggrega- 
tion, in which the particles trajectories are in a single 
direction. They conclude that the 7 exponent is nonuni- 
versal, i.e., depends on the lattice structure. 

In Fig. the double-logarithm plots of p — p as a 
function of r for distinct A values are shown. The den- 
sity po was obtained by searching for the best linear fit 




FIG. 5: Double-logarithm plots of p — po against r for A = 0.3 
(squares) and A = 0.99 (circles). The dashed line corresponds 
to the slope —0.46. The linear fits of the data provide 7 ~ 0.45 
for A = 0.3, and 7 « 0.47 for A = 0.99. 



in the larger linear region. To avoid the active region, 
we limited the fits to those data corresponding to a half 
of the cluster sizes. Depending on the A value, lattices 
with linear size L = 5 x I0 3 or L = 10 4 , and 10 to 20 
independent runs were used. One can observe a power 
law regime for r > 10 showing that the approach to the 
stationary value obeys Eq. [5] In Fig. the asymptotic 
density po and the 7 exponent are shown as a function of 
1 — A, the distance from the transition point, po acts as 
an order parameter, which vanishes at the critical point 
following a relation po ~ |1 — A|^. The exponent obtained 
from the data of Fig. E{a) was /3 = 0.27(1), whereas the 
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exponents obtained for L = 2 x 10 3 and L — 10 4 were 
13 = 0.28(2) and (3 = 0.26(1), respectively. The num- 
bers in parenthesis represent the uncertainties. In Fig. 
E^b), the 7 exponents for distinct A values are shown. 
One can observe that 7 fluctuates around the value 0.46. 
The smaller value found was 7 ps 0.43 and the larger one 
7 w 0.48. Our simulations suggests that the 7 exponent 
is independent of A and its value is 7 = 0.46 ± 0.02. The 
error indicated in the 7 value was evaluated through an 
average over the data of Fig. HJb) . 

In order to test the universality of the 7 and (3 expo- 
nents we studied two versions of the present model. In 
the first one, we used square lattices, but the walkers stick 
to the cluster if they reach a nearest or a next-nearest 
empty neighbor of an occupied site. In the second one, 
we used a hexagonal lattice. Lattices with size L = 2000 
were used. The exponents for the first modified version 
were = 0.27 ± 0.02 and 7 = 0.49 ± 0.03, and the ex- 
ponents for hexagonal lattice were (3 — 0.28 ± 0.02 and 
7 = 0.49 ± 0.03. These results lead us to conclude that 
these exponents are universal. 




10~ 2 10 1 10° -m 2 10 1 10° 

FIG. 6: (a) Stationary density po and (b) 7 exponent denned 
in Eq. ® as a function of the distance from the transition 
point. The lattice size used was L = 5 x 10 3 and simulation 
stopped when the cluster reaches the lattice edge. 

The morphological transition between BA and DLA 
models was also characterized by the crossover radius £ 
defined in Eqs. Q and @. The number of particles 
M(r) inside a region delimited by a circle of radius r 
centered at origin was evaluated. The M vs. r curves 
exhibit tenuous crossovers determining the transition be- 
tween DLA- and BA-like scaling regimes. In Fig. EJa), an 
example of this crossover is shown. Since the growth pat- 
terns scale as DLA (BA) for small (large) length scales, 
in order to evaluate the crossover £ we fitted the curves 
by power laws M(r) ~ r df , where df — 1.71 and df = 2 
were used for the initial and the final curve regions, re- 
spectively. The crossover lengths obtained through this 
method are drawn as a function of the distance from 
the transition point in Fig. [7Jb). The length £ diverges 



at A = 1 following a power law £ ~ |1 — A| _I/ , where 
v = 0.61(1). Moreover, the mass at the critical point 
diverges as ~ |1 — A|~ Q , where a = 0.97(2). 

As discussed in section[|] these crossovers between frac- 
tal and homogeneous patterns occur due the crossover in 
the particle trajectories. However, the crossover length 
of the walker trajectories is given by (see appendix) 

ttA sin(TrA) 

t,w = . -, rr : — , (6) 

sm(7rA) 7rA 

which diverges as £ w ~ |1 — A| _1 for A < 1. Thus, al- 
though the transition between DLA and BA models is 
due to the transition in the walk dimensionality, the re- 
spective crossover lengths are not proportional. 




r 



FIG. 7: (a) Determination of the crossover between DLA- 
and BA-like scaling regimes for A = 0.90. The straight lines 
represent the slopes 1.71 and 2, i.e., the fractal dimensions for 
DLA and BA models, respectively, (b) The crossover length 
(squares) and the correspondent mass (circles) as a function 
of the distance from the transition point. These results were 
obtained for L = 5 x 10 3 . 

Notice that Eq. JHJ describes the mean density behav- 
ior when r <J £. Moreover, p ~ r df ~ 2 when r <§; £ due 
to the cluster fractality in this length scale. Thus, using 
Eq. JSJ), we found that the mean density at the crossover 
can be written as 

^=yli|l-A|^ + i4a|l-Ar, (7) 

where A\ and A2 are constants. But, the mean density 
at the crossover is given by 

p^Mi~\i-\\^-«. (8 ) 

Comparing Eqs. and ©, we have that they are con- 
sistent only if 

(3 = isj and (3 = 2v - a. (9) 

In agreement with the scaling relation ©, the number of 
independent exponents are reduced from 4 to 2. Using 
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the exponents measured for systems with size L = 5000, 
we found v-f = 0.28(2) and 2v - a = 0.25(3), beside 
(3 = 0.27(1). The difference between these values is inside 
of the error margins indicated in the parenthesis. The 
large uncertainities obtained in the exponents (5%-10%) 
are originated in the difficulty in the determination of the 
exact crossover points. 



IV. CONCLUSIONS 

In the present work, we studied the transition between 
diffusion-limited aggregation (DLA) and ballistic aggre- 
gation (BA) models. We used a model in which the ran- 
dom walks in the DLA model are replaced by biased ran- 
dom walks with a drift in a random direction. The drift is 
controlled by a parameter A 6 [0, 1] that leads the model 
from BA (A = 0) to DLA (A = 1) (see Eq. ©). Also, an 
efficient algorithm, which allows large scaling analysis of 
the growth patterns, was introduced. 

For any bias, the clusters are fractal (DLA-like) on the 
short length scales whereas nonfractal patterns are ob- 
tained on the large ones. The transition between DLA- 
and BA-like scaling regimes is determined by a char- 
acteristic length £ that diverges as A — > 1 following a 
power law £ ~ |1 — A| _l/ , where v = 0.61(l),(r) while 
the cluster mass at the crossover follows the relation 
M;; ~ |l — X\~ a , where a — 0.97(2). This crossover was 
not numerically determined in similar previous works. 
The density in the inner regions of the cluster reaches 
an asymptotic value po ^ 1 1 — A|^, where (3 — 0.26(1). 
However, this approach is slow and follows a power law 
decay with a universal exponent 7 = 0.46(2) independent 
on the drift. These exponents obey the scaling relations 
/3 = 1/7 and = 2v — a. 

It is worth to stress two main contributions of the 
present work. The first one is the development of an 
algorithm that can be used to study other models with 
biased random walks as for example those related to de- 
position processes |l5l Il6[ , for which the determination 
of universality classes is a hard work. The second one 
is the careful quantitative characterization of the transi- 
tion between DLA and BA growth models that, in our 
knowledge, was not previously done. The understanding 
of these crossovers can be an essential tool in the analysis 
of real fractals, which always exhibit scaling on limited 
ranges. 
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APPENDIX: DEMONSTRATION OF EQ. © 

For sake of simplicity, we consider Eq. Q with a drift 
direction tp = and xq = yg = 0. Iterating Eq. J2J) for n 
steps, we found 



(x n ) = y^(cos(Agj)), 



i=l 



and 



»=i j=i 



But, 



(cos(A0<)) = — / cos(X9)d9 = sm ( nX > , 

27T J_„ 7rA 



(A.l) 



(A.2) 



(A.3) 



and 



(cos(A#i) cos(A#j)} 



sin(7rA) 

ttA 
sin(27rA) 



2ttA 



(1-%) 



+ 1 



(A.4) 



where Sij is the Kronecker delta. 

Substituting Eqs. (|A.3jl and (|A.4jl in Eqs. IjA.lj) and 
l2j, respectively, we found 



(x n ) = n 



sin(7rA) 
ttA 



and 



a 2 x (n) = n- 



sin(27rA) 




sin(7rA) 




2ttA 




ttA 





(A.5) 



(A.6) 



where cr%(ri) = (a;„) — (x n ) 2 is the variance of the coor- 
dinate x n . With a similar analysis, we obtained 



and 



(Vn) = 

sin(27rA) 



2ttA 



(A.7) 



(A. 



Thus, the walk mean displacement r(n) = 
V ( x n) 2 + (Un) 2 and the variance a 2 (n) 
are given by 



o2 + vl 



r(n) = n 



sin(7rA) 
ttA 



and 



<r(n) = n 



sin(7rA) 



TTA 



1/2 



(A.9) 



(A.10) 
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The crossover of the walk dimensionality (d = 2 for 
short times and d = 1 for long ones) occurs when r(n) ~ 
<r(n). Making equal Eqs. 1|A.9J| and (|A.10|I . we obtained 
an estimative of the characteristic number of steps Af 
necessary for the crossover 



FIG. 8: Double- logarithm plot of the crossover length £ w as a 
function of 1 — A. The solid line represents the data obtained 
from Eq. HA. 121 and the dashed line has a slope —1. 



N ■ 



1. 



sin(7rA) 

Thus, the characteristic crossover length is 



-(AO 



7rA sin(7rA) 



sin(7rA) 7rA 



(A.ll) 



(A.12) 



In Fig. ijSJ), £ w is plotted as a function of 1— A. Expanding 
Eq. QA. 12(1 around A = 1, we found £ w ~ |1 - A| _1 . 
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